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ABSTRACT 

We present deep GMOS-S spectroscopy for 11 galaxy groups at 0.8 < z < 1.0, for galax- 
ies with tab < 24.75. Our sample is highly complete (> 66%) for eight of the eleven 
groups. Using an optical-NIR colour-colour diagram, the galaxies in the sample were sep- 
arated with a dust insensitive method into three categories: passive (red), star-forming (blue), 
and intermediate (green). The strongest environmental dependence is observed in the frac- 
tion of passive galaxies, which make up only ~ 20 per cent of the field in the mass range 
10 10 3 < Mstar/Af© < 10 110 , but are the dominant component of groups. If we assume 
that the properties of the field are similar to those of the 'pre-accreted' population, the envi- 
ronment quenching efficiency (e p ) is defined as the fraction of field galaxies required to be 
quenched in order to match the observed red fraction inside groups. The efficiency obtained 
is ~ 0.4, similar to its value in intermediate-density environments locally. While green (in- 
termediate) galaxies represent ~ 20 per cent of the star-forming population in both the group 
and field, at all stellar masses, the average sSFR of the group population is lower by a factor 
of ~ 3. The green population does not show strong H<5 absorption that is characteristic of star- 
burst galaxies. Finally, the high fraction of passive galaxies in groups, when combined with 
satellite accretion models, require that most accreted galaxies have been affected by their en- 
vironment. Thus, any delay between accretion and the onset of truncation of star formation (r) 
must be < 2 Gyr, shorter than the 3 — 7 Gyr required to fit data at z = 0. The relatively small 
fraction of intermediate galaxies requires that the actual quenching process occurs quickly, 
with an exponential decay timescale of r q < 1 Gyr. 

Key words: galaxies: evolution - galaxies: luminosity function, mass function - galaxies: 
general 



1 INTRODUCTION 



Studies of galaxy evolution have shown that the star formatio n 
density of the uni verse depends on epoch (e.g. iLillv et al.|[l996h . 
stellar mass (e.g. iBrinchmann et al.| 20041) . and environment (e.g. 



iBaldrv et al] 120061 : ISobral et alj 1201 ll) . At z ~ 1, star-forming 



galaxies continue to show a strong relationship between the specific 
star formation rate ( sSFR) and stellar mass, with moderate scatter 
dNoeske et alj2007l : ISalmi et alj|2012h . Galaxies that lie above this 



sequence tend to either host an active galactic nucleus (AGN) or a 
short-lived burst of star formation, while those below the relation 
may be cases where th e galaxy's star formation is shutting down 
dWhitaker et al.l2012bh . 

Groups are an important test of the effects of environment 
on gala xy evolution, as the majority of galaxies are found inside 
groups ( lEkeetalj|2004h and they are a crucial link in the hier- 
archical growth of structures from individual galaxies to dense 
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clusters. Galaxy groups have been studied i n the context of gen- 
eral redshift surveys, for example . CNOC2 Jcarl berg et aL bOOll : 
iFinoguenov e7ai]|2009h. D EEP2 dGerke et alj|2012h . and zCOS- 



MOS (Knobeletal. 2012a), including imaging- only photometric 
z wo rk (e.g. iGiodini etal] 1201 2| ; iLi et al.l 120121; iLeauthaud et alj 



2012h . Another approach is to study individual grou ps in more 
detail and with highly complete spectrosc opy (e.g. Ijeltema et alj 
120091 IPello etai]|2009l ; lBalogh et alj201lh . 

Observations have also demonstrated that there already exists 
a strong correlation between gal axy environment and the fraction 
of pa ssive galaxies at z ~ 1 (e.g. lCooper et al.l20 06; Pozzet tiet alj 
l20ld) . In general, the fraction of quiescent galaxies is always 
highe r in denser and more massive systems dvan den B osch et al. 
20081: iGeorge et alj|201 \Y iMuzzin et al]|2012l: iTanaka et alj|2012l; 



Presotto et al.ll2012h . ex cept where the galaxy merger rate ma y be 
exceptionally high (e.g. iGeach et al]|201ll ; ISobral et alll201 lh . At 
2 > 0.8, there is evidence that, at least in some cases, the SFR 
of the star-forming population in l ocally-dense environments may 
in fact be higher than in the field dCooper et al.l l2008 : Elba z et al] 
20071 : ISobral et alj|201 H ; lli et al.ll201 ll). while o ther groups have 
not found this enhancement (e.g. lPatel et alj[2oTTh . 

It is perhaps most natural to interpret these observations in 
the context of the halo model, in which the properties of the 
host dark matter halo completely determine the galaxy's proper- 
ties and where the evolution of satellite galaxies can be distin- 
guished from that of central galaxies. In this context, the radial 
gradient s of the quiescent f r action observed in gr o ups and clus- 
ter (e .g. iBalogh et all l2000l: lEllingson etai]|200ll; IPresotto et all 
2012h reflect the sate llites' accretion history (e.g. lSmith et alj2012l ; 



De L ucia et alj2012h . 



is inefficient. By trying to match the full shap e of the observed 
SFR distribution at z = 0, Wet zel et al.1 J2012h recently proposed 



An empirical model of galaxy evolution that accounts for most 
of these trends was put forward bv lPeng et al.|j2010ll2012l) . In this 
model, galaxy evolution is driven by two "quenching rates", one 
that depends, perhaps indirectly, on stellar mass and another that 
is related to environment (local density). While the efficiency of 
environmental quenching in this model is roughly independent of 
stellar mass and redshift, its effects are most prominent amongst 
low-mass galaxies, M < 1O 1O 5 M , for which mass-quenching 

hape 
I3)r 

that this environmental dependence must take the form of a rapid 
(r < 1 Gyr) satellite quenching timescale, but preceded by a long 
delay of 2 — 4 Gyr. 

Both of these models make predictions about the environmen- 
tal dependence of galaxies at z = 1, which are poorly constrained 
by existing observations, in part because the effects are expected 
to be most visible in the low-mass star-forming population. In par- 
ticular, such a long delay before a satellite begins to respond to its 
changed environment means that trends with environment should 
be much weaker at z ~ 1, since at that epoch m ost galaxies 
would not have spent a long enough time as a satellite dMcGee et alj 
l2009h . 

For all of the reasons listed above, the study of the universe at 
0.8 < z < 1 is particularly important. Compared to z ~ 0, this 
epoch has high in-fall rates into groups, and the galaxies themselves 
have high star formation rates and gas fractions, thus p resenting the 
best o pportunity to observe any transitional objects (McG eeet alj 
l2009h . 

This paper will present the data collected as part of the 
Group Environmen t Evolution Collaboration 2 (GEEC2) study 
(Balo gh~et al.l20 1 lh . Our survey selects a sample of ~ 20 robustly- 
identified groups with extended X-ray emissions (Finogueno v et alj 



120071: IGeorge et al| l201lh from the COSMOS survey fields 
dScoville et alJl2007T> . Then, we seek to obtain deep, highly com- 
plete spectroscopy for their member galaxies, focusing on low 
mass galaxies down to 10 9 Mq, where environmental effects may 
be most prominent. This approach of surveying a smaller num- 
ber of well-sampled groups h as been previously applied at z 



ippji 

0.5 with great success (e.g. IWilman et al" 120051: IBalogh et al] 



20071: iMcGee et al]|2008l; 



Wilman ct al 



Finog uenov et alj2009l : IConnellv et alj|2012h 



20091 : IBalogh et al]|2009t 



Analysis of low-mass groups at z ~ 1.0 has reve aled a distinct 
popul ation of intermediate-colour (green) galaxies (Balo gh et al] 
1201 lh . which may indicate the presence of a transition popula- 
tion. These green galaxies have morphologies half-way between 
the blue-cloud and red-sequence galaxies and do not appear to be 
exceptionally dusty. Inside the group environment, they may be ev- 
idence for the impact of environment on star formation. The sample 
has since been expanded with additional GMOS observations, al- 
lowing us to improve the statistical significance of these trends, and 
study their SFR and stellar mass properties. 

The structure of the paper will be as follows. In § 2, we present 
a discussion of the updated GEEC2 sample with new observations 
from the 11A semester and the updated data reduction process. 
In § 3, the main analysis of the data will be presented, includ- 
ing our separation of the sample into sub-populations using their 
colour, the determination of completeness weightings, and the cal- 
culation of star-formation rates. In § 4, the results are presented, 
including stellar mass functions and the comparison of sSFR and 
H<5 strengths for different sub-populations and environment. Fi- 
nally, § 5 will include a discussion of how satellite accretion models 
and the red/green/blue fractions in our sample constrain quenching 
models and the effects of environment on the group population. A 
companion paper (Hou et al., in prep) considers how the galaxy 
population in these groups depend on the group dynamics. 

Throughout the paper, a cosmology with Q m = 0.272, Aa = 
1 — n m and h = H a / (lOOkm/s/Mpc) = 0.702 was assumed 
dKomatsuetai]|201ll) . All magnitudes are in the AB system and 
star formation rates are given for the Chabrier IMF. 



2 OBSERVATIONS 

A more complete discussion a nd explanation of the GEEC2 survey 
can be found in IBalogh etalj d201ll) . which presents results from 
the first semester of Gemini data. The survey is designed to cover 
20 groups at 0.8 < z < 1, requiring ~ 160 hours of Gemini time 
over three years. So far, a total of 82.5 hours of data have been col- 
lected from Semesters 10A and 11 A, covering approximately 11 
groups. The following is a summary of the most important details, 
including an update of the sample to include observations from 
Semester 1 1A. 



2.1 Galaxy Group Selection 

The GEEC2 su rvey targets galaxy groups identified within the 
COSMOS field ( Scoville et al. 2003 ) based on an alysis of the deep 
X-ray data dFinoguenov et al.ll2007l. [20091 FzOlO). Identification of 
group redshifts is done using all available spectroscopy, including 

but not limited to, the 10K zCOSMOS survey 

W e considered all groups presented in the lGeorge et al.l d201 lL 
l2012h catalogue that lie within the redshift range 0.8 < z < 
1, have at least three spectroscopically-determined members, and 
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were classified as either category 1 (good detection, with a well- 
determined X-ray centre) or category 2 (secure detection, but with 
an unreliable X-ray centre). There are 21 groups satisfying this se- 
lection. Of these, we have given preference to the lowest-mass, 
highest redshift systems. The X-ray selected group s in this paper 
are a subset of the groups in the lKnobel etai1d2012d) spectroscopic 
catalogue. 

Over two semesters with Gemini we observed ten of these 
groups, and from these data discovered an additional serendipitous 
group in the same redshift range (213a). Their properties are tabu- 
lated in Tabled 



2.2 Gemini Observations 

2.2.1 Spectroscopic Target Selection 

Targets were sele cted from the COSMOS photometric catalogue 
(Capak et alj|2007l) . All optical colours and magnitudes are com- 
puted within a 3" diameter aperture. A crucial part of our strategy 
is the use of photometric redshifts to select potential group mem- 
bers. Precise photometric redshifts are derived from th e photome- 
try of 30 broad, intermediate and narrow-band filters (Ilbe rt et alj 
2009), using a template-fitting technique, and calibrated based on 
large spectroscopic samples. Even for the faintest galaxies in our 
sample, most of the objects have photometric redshifts determined 
to better than 10 per cent. 

We give highest priority to galaxies that have 21.5 < r < 
24.75, and a redshift within 2oz P hot of the estimated group red- 
shift, where a zp hot is the 68 per cent confidence level on the pho- 
tometric redshiftQ Secondary priority slits are allocated to galaxies 
with 15 < r < 24.75 and 0.7 < z pho t < 1.5. 

We designed 3 — 4 GMOS masks on each target, using the 
GMMPS software. On each mask we are able to allocate 40 — 50 
slits. Typically, well over half the slits on the first mask are allocated 
to our top priority objects. This fraction decreases on subsequent 
masks, as we use up these targets. In most cases, with three masks 
we are able to target at least 40 of these high priority galaxies. Many 
masks for a given target use some of the same alignment stars and 
are always at the same position angle; thus a small fraction of the 
CCD area is unusable regardless of the number of masks obtained. 



2.2.2 11A Gemini Observations 

Our 10A observations were described in Balogh et al. I d201ll) . In 
11A we obtained an additional 32.5 hours of data from Gemini 
South, in the Band 1 queue. This returned science data of 2 hours 
on-source exposure, for 12 masks in 4 groups. All science obser- 
vations were obtained in clear conditions with seeing 0.8 arcsec or 
better in i. 

The sp ec trosco pic setup was almost the same as in 
IBalogh et all d201 lh. We used th e nod & shuffle mode 
dGlazebrook & Bland-Hawthornl 2001), observing the galaxy in 
both nodded positions. Slits were 1 arcsec wide, and we use the 
R600 grism with the OG515 order blocking filter. One difference 
from the 10A observations is the detector was binned by 4 in the 
spectral direction, rather than 2. This yielded a dispersion of 1.86 



1 At our redshift of interest, the average uncertainty on Zphot increases 
from cr zvhot ~ 0.007 for the brightest galaxies to CT zp ^ ot ~ 0.04 for 
those at our limit of r = 24.75. Even for these faintest objects, 90 per cent 
of the galaxies have z p hot uncertainties of less than <T zphot < 0.07. 



A/pix, and a spectral resolution, limited by the slit width, of ~ 12.; 



2.2.3 Data reduction 

All data are reduced in IRAF, usin g the GEMINI packag es with mi- 
nor modifications as described in iBalogh et all d201ll Redshifts 
are measured by adapting the ZSPEC softwa re, kindly provided b y 
R. Yan, used by the DEEP2 redshift survey dNewman et alj|2012h . 
This performs a cross-correlation on the ID extracted spectra, using 
linear combinations of template spectra. The corresponding vari- 
ance vectors are used to weight the cross-correlation. Finally, red- 
shifts are adjusted to the local standard of rest using the IRAF task 
RVCORRECT, though this correction is negligible. 

We adopt a simple, four-class method to quantify our redshift 
quality. Quality class 4 is assigned to galaxies with certain red- 
shifts. Generally this is reserved for galaxies with multiple, robust 
features. Quality class 3 are also very reliable redshifts, and we 
expect most of them to be correct. These include galaxies with a 
good match to Ca H&K for example, but no obvious corroborating 
feature. 

Note our spectral resolution was a factor of two lower in our 
1 1 A observations, and this meant that it was generally not possible 
to resolve the [Oil] doublet. Thus, single emission lines are gen- 
erally assumed to be [Oil] and given quality class 3. This is, in 
part, due to our prior photo-z data for these objects. We would also 
use this in cases where H&K are detected in a region of telluric 
absorption, but there is at least one other likely match to an absorp- 
tion feature. We take particular care not to assign class 3 or 4 to a 
galaxy for which H&K are the only identifiable features, and lie on 
a telluric absorption line. 

Class 2 objects correspond to "possible" redshifts. These in- 
clude some spectra that are reasonably likely to be correct (e.g. 
H&K on top of a telluric line and no other corroborating features), 
but also some that are little more than guesses. Class 1 is reserved 
for "junk", with no chance of obtaining a redshift. In this analysis 
we only consider galaxies with class 3 or 4 quality redshifts, which 
represent 603 of the 810 unique galaxies for which a good spectrum 
was extracted. 

We check our redshift accuracy in two ways. There are 
ten galaxies with redundant observations, where both observation 
yielded a class 3 or 4 redshift. The dispersion of the redshift dif- 
ferences for these objects is ~ 150 km/s, meaning the typical 
redshift uncertainty is ~ 100 km/s. Secondly, we deliberately re- 
o bserved 1 1 galaxies w ith redshifts available from zCOSMOS . As 
in IBalogh et all feoill) . we find a small systematic offset, such that 
our redshifts are smaller than those in zCOSMOS. Over the two 
semesters we have 19 duplicates of this type. Omitting one outlier, 
the bias is Az = 6.2 x 10~ 4 , corresponding to a rest-frame veloc- 
ity bias of ~ 100 km/s at z = 0.9. Then, we adjust the zCOSMOS 
redshifts by this small amount so they are consistent with ours. 



2.2.4 zCOSMOS 

We also make use of the DR2 release of t he 10K zCOSMOS spec- 
troscopic survey dLillv et all I2007L l2009h . from which we obtain 
spectra, photometry and redshifts for moderately bright galaxies 
(i < 22.5) over most of the survey area. This is a sparsely-sampled 
redshift survey, with ~ 40 per cent sampling completeness. All 
galaxies with redshift quality greater than 2.0, and with a high prob- 
ability (> 90 per cent) of being correct (note the zCOSMOS quality 
flags are defined differently from our own) are used. 
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Group 


RA 


Dec 




N 
* mask 


N 


(j (km/s) 


R, mo (Mdc) 


Comp. (%) 


M,i„„(10 13 Mr>) 


40 


150.414 


1.848 


0.9713 


2 


15 


690 ±110 


0.34 ±0.04 


64 


11 ±4.7 


71 


150.369 


1.999 


0.8277 


3 


21 


360 ± 40 


0.34 ±0.03 


82 


3.0 ±0.9 


120 


150.505 


2.225 


0.8358 


3 


3 1 


/ton _1_ ar\ 
480 ± 60 


n *7Q _l_ n n*7 
0. (a ± 0.0/ 


to 
OV 


LZ. ( ± 4.2 


121 


150.161 


2.137 


0.8373 


3 


5 


210 ± 60 


0.09 ±0.01 


100 


0.2 ±0.2 


130 


150.024 


2.203 


0.9374 


3 


34 


600 ± 70 


0.71 ±0.06 


76 


17.8 ±5.7 


134 


149.65 


2.209 


0.9467 


3 


23 


450 ± 60 


0.97 ±0.07 


70 


13.3 ±4.7 


143 


150.215 


2.28 


0.881 


3 


20 


580 ± 60 


0.23 ±0.03 


100 


5.2 ± 1.7 


150 


149.983 


2.317 


0.9334 


4 


25 


300 ± 40 


0.89 ± 0.08 


75 


5.5 ± 2 


161 


149.953 


2.342 


0.944 


4 


8 


170 ± 30 


0.53 ±0.12 


70 


1.0 ±0.5 


213 


150.41 


2.512 


0.879 


2 


9 


260 ± 100 


0.84 ±0.13 


52 


3.8 ±3.4 


213a 


150.428 


2.505 


0.9256 


2 


8 


110 ± 30 


0.62 ±0.09 


40 


0.4 ±0.3 



Table 1. Properties of the eleven galaxy groups observed with GMOS in semesters 10A and 1 1A. The position, median redshift z me d> 
rest-frame velocity dispersion, and the number of group members are determined from our GMOS spectroscopy (combined with 
available zCOSMOS 10k data), as described in the text. The number of GMOS masks observed in each field is given by -/V ma sk; note 
that groups 150 and 161 lie within the same field, as do groups 213 and 213a (a serendipitous discovery in the background). iJrms 
is the rms projected distance of all group members from the centre. The penultimate column lists the spectroscopic completeness of 
the group, where we list the percentage of candidates within i? r ms of the group centre for which spectroscopic redshift was obtained, 
where a candidate is a galaxy with a photometric redshift consistent with the group redshift at the 2<r level (see § 13.4.11 . Finally we 
list the dynamical mass of the group, calculated from a and -R rms following Equation[T] 

3 ANALYSIS 



3.1 Group Masses and Membership 

The groups in our sample have already been robustly identified 
from sparse spectroscopy and deep X-ray observations. With the 
additional spectroscopy we do not therefore start with a group- 
finding algorithm, but only improve existing estimates of the cen- 
troid and the velocity dispersion. As these are confirmed, sig- 
nificant over-densities with extended X-ray emission and well- 
determined centres, the amount of possible field contamination 
should be substantially lower than the ~ 15 per cent exp ected from 
a general group-finding algorithm (Knob eletal1l2012al) . 

For each group, we start by considering all galaxies within 
r — 1 Mpc and 4000 km/s (rest-frame) of the nominal X-ray group 
centre from lGeorge et al.l d201 ll 12012). The velocity d ispersion of 
these galaxies is determined using the gapper estimate (Be ers etlil] 
1990). We also compute the (unweighted) mean spatial position 
of the galaxies, and the rms projected separation from this centre, 
which we call R rms - We then iterate this process, clipping galaxies 
with velocities v > C z a from the median redshift, and positions 
r > CrRrms from the recomputed centre. 

For most groups, we use C z — 2 and CV = 2. Groups 150 
and 161 are close together in redshift and position, and we adjusted 
these parameters to better separate them. We adopt C z = 1.5 and 
C r = 2 for group 150, and C z = 1.4 and CV = 1.6 for group 161. 
Finally, groups 121 and 213a have few members and we needed to 
relax CV. = 2.5. 

We adopt the final a and i? rms from this procedure, and de- 
fine group members as those within 2.5<r of the median redshift 
and within r < 2R rm3 , although the limiting radius can be varied 
depending on the application. Uncertainties on a and i? rms are then 
computed using a jackknife method, iterating only over this list of 
group members. Thus, these uncertainties do not include system- 
atic biases due to the c lipping process, wh ich are likely to be at 
least ~ 15 per cent (e.g. lWilman et alj|2005l) . 

The final velocity dispersions, reported in Table [T] are cor- 
rected to their rest-frame values. No explicit correction is made for 
the effect of redshift uncertainties on these dispersions. The disper- 
sions of the few groups with a < 200 km/s are thus likely biased 




Figure 1. For confirmed group members with 0.7 < z p hot < 1-5 we 
show the difference between spectroscopic and photometric redshifts, as a 
function of (V — z) colour. Asymmetric error bars are 68 per cent confi- 
dence intervals on Zphot ■ Filled symbols are quality class 4, which means 
the spectroscopic redshift is certain, while open symbols are class 3 and 
thus generally reliable (and used throughout this analysis). The red points 
identify those few group galaxies that were allocated lower priority during 
mask design because their 2 p hot i s more than 2<r away from the mean group 
redshift. These make a negligible contribution to the final membership. 



For most of the analysis, we only consider galaxies within the 
field-of-view of our GMOS fields, which simplifies and homoge- 
nizes the selection function. In some cases where this is not impor- 
tant, however, we consider the whole catalogue. 
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high, perhaps up to a factor of ~ 2 dKnobel et ai] |2009). Note that 
this systematic bias can be larger than the statistical uncertainties 
on the velocity dispersion. 

Group dynamica l masses are computed from a and i? rms , as 
in lBdoghetailteOllh : 



Mdyn — firms d 



(1) 



These masses, listed in Table Q] range from ~2x 10 12 M© to 
~ 1.8 x 1O 14 M . 

Any galaxy at 0.8 < z < 1, within the GMOS fields-of-view, 
but not assigned to a group, forms part of our comparison sam- 
ple. This is not necessarily a low-density sample, as it may contain 
unidentified groups, but nor is it representative of the overall field, 
since some groups have been removed. For this paper, we will usu- 
ally refer to it as the "field" sample for simplicitjQ Any difference 
from a truly unbiased field population is probably small, at least 
compared with our group sample. 



3.2 Stellar Mass Measurements and k-corrections 

We fit the spectral energy distribution for each galaxy us- 
ing all available ph otometry, following the meth od descr i bed in 
iMcGee et all < l201ll) . but now using the updated iBruzuall fc007l) 
models, commonly referred to as Chariot & Bruzu al (2007) or 
CB07. As previously, we assume a lChabrien (2003) initial mass 
functio n and, following the general procedure of [Salim et alj 
d2007h . assume the star formation history of a galaxy can be repre- 
sented by an exponential decay model with superposed bursts. 

Every model in the template suite is compared with the data, 
and the \ 2 difference is computed. This is then used to define a 



weight uii = exp (-x?/2 



for model i, and we construct a prob- 
ability distribution function (PDF) for every parameter of interest 
by compounding these weights. We adopt the median of the PDF 
as our best estimate of the parameter, and the 16 — 84 percentile 
range of the PDF is taken to represent the ±lcr uncertainty. 

We k-correct all colours to z = 0.9, the redshift of in- 
terest for our survey, u sing the KCORRECT IDL software of 
iBlanton & Roweisl fc007t) . to minimize the effect of such correc- 
tions. We denote these colours as, for example, (V — z) 0,9 . The 



perfectly straddles the 4000 A break at 0.8 < z < 1. 



It is roughly comparable to rest-frame (u — B), but in fact, the V 
filter is measuring rest-frame flux intermediate between GALEX 
NUV and ground-based u. 

The difference betw een the stellar masses based on the older 
iBruzual & Charlotl d2003l) models, and the new ones, are generally 
small Jllbert et aljioich . However, we have made some other, more 
important chang es to how stellar ma sses are computed since our 
previous work in Bal ogh et al.1 d201ll) . First, we apply a correction 
to the aperture magnitudes to approximate total magnit udes, which 
leads t o a systematic increase in stellar mass. Following lllbert et alj 
(2010), we also test the sensitivity of the \ 2 distribution to small 
zero-point offsets, and the omission of particular filters. We find 
that removing the J-band data (where strong Ha emission is often 
present), and applying small zero-point offsets, significantly im- 
proves the number of \ 2 outliers, and thus we adopt this modifica- 
tion to the fits. 



2 In lBaloghetailfcOllI) we called this the "non-Xgroup" sample but we 
decided to discard this terminology here in favour of the "field" sample. 
There is no change in its definition. 



Index 


Blue Continuum 


Line Definition 


Red Continuum 


[On] 


3653 - 3713 


3722 - 3733 


3741 - 3801 


H<5 


4030 - 4082 


4082 - 4122 


4122-4170 



Table 2. The table list the wavelength intervals used for the estimation of the 
blue and red continuum, as well as the line definition for the spectroscopic 
indices used in this paper. The values in the table are provided in units of 
Angstroms. 



This procedure leads to masses that are ~ 50 per cent larger on 
average from those we used in Balogh et al 1), with a standard 
deviation of about a factor ~ 2. 



3.3 Spectroscopic Indices 

Spectral line indices were measured for all galaxies with secure 
redshifts. First, the ID spectra are shifted into each galaxy's rest- 
frame with the IRAF DOPCOR command. The FITS files are then 
imported into IDL, where the relevant spectroscopic indices can 
be measured using a custom-written script. The equivalent widths 
of the [On] and H<5 lines are defined by first fitting a straight line 
between the centres of the red and blue continuum ranges, where 
the estimates for flux levels of the red and blue continuum were ob- 
tained by calculating the average flux over the wavelength ranges 
presented in Table [2] Next, the flux over the line index is summed 
over the wavelength ranges listed in Table |2]for each spectral fea- 
ture. Using the calculated linear fit to the continuum, the equivalent 
width for the H<5 and [On] lines can be determined using the equa- 
tion: 



Wa = 



1-| MA. 



(2) 



where F c is the flux level of the continuum and Fa is the total flux 
per unit wavelength. 

For this paper, we adopt the sign convention such that a pos- 
itive [Oil] equivalent width indicates an emission feature, while a 
positive H<5 equivalent width corresponds to an absorption feature. 



3.3.1 Weighted-Median and Stacked-Spectra Methods 

Since the signal to noise ratio for these faint spectra are relatively 
low (S/N ~ 1 per pixel), any measurements of spectroscopic 
indices (including their derived quantities, like the star formation 
rate) can have considerable uncertainty. It is therefore valuable to 
consider quantities averaged over various subpopulations, and we 
do this in two different ways in this paper. 

For the first method, the relevant indices for each galaxy spec- 
trum are measured individually, using the procedure outlined in 
§ 13,61 We then compute a median that includes the galaxies' spec- 
troscopic and mass-completeness weights (the derivation of these 
weights will be further described in § I3.4.U . For this weighted- 
median method, the relevant quantities (such as [Oil] luminosity) 
are first sorted in order of increasing value. Then, the total weight is 
calculated. The point at which the cumulative weighting is equal to 
half of the total weight is considered to be the weighted-median 
value. The uncertainties in this method can be determined with 
the bootstrap method. The distribution is randomly re-sampled 500 
times and the median value re-calculated for each iteration. The 
standard deviation of this distribution of median values can provide 
a robust error estimate for the original measurement. 
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The second method is to stack the individual spectra and mea- 
sure the feature on the combined spectrum. First, the spectra are 
multiplied by their total weights, which is a combination of their 
spectroscopic and mass-completeness weighting. Then, the spectra 
are combined by taking the average flux of each pixel, using the 
IRAF function SCOMBINE. The resulting stacked spectra is used 
for the subsequent calculations for the properties of this bin. The 
measurement errors in the equivalent width s of the stacked spectra 
can be calculated using the formula from Vollmann & Eversberd 
J2006h . 

a(W A ) = yi + ¥ ^ 7 ^ (3) 

where the S/N term is the total signal to noise ratio over the line 
index in question. This is defined as the mean value of the flux per 
pixel, divided by the variance and multiplied by the square root 
of the number of pixels in the line index definition. Note that this 
estimate represents only the uncertainty in the mean value, and does 
not include the variance associated with the distribution of pixel 
values entering the stacking process. 

3.3.2 Active Galactic Nuclei 

We identify strong AGN by cross correlating with the point source 
catalogue from deep Chandra observations within 5". Only seven 
galaxies have such a match, and they are identified in all appro- 
priate figures as symbols with a thick, black square around them. 
Weak, narrow-line AGN, or heavily-obscured AGN could still be 
present in the sample. Radio-loud AGN could also be present, but 
these are not expected to bias SFR measurements. 




(V-z) 09 

Figure 2. We show the stellar mass-to-light ratio as a function of colour, 
for galaxies with secure redshifts in our sample. The mass is derived from 
SED fitting using lBruzua] 120071) models with a lChabrierl d2003l) IMF, and 
the luminosity here does not include a k-correction. Black points represent 
galaxies at the lower end of our redshift range, 0.8 < z < 0.85, while the 
red points correspond to the upper end at 0.95 < z < 1.0. The solid and 
dotted lines (barely distinguishable over most of the range) are fits to the 
median M/L r for the low and high redshift samples, respectively. The fit 
is bilinear, with a different slope for red and blue galaxies. 



3.4 Completeness Corrections 

3.4.1 Spectroscopic and Redshift Completeness 

There are several potential sources of incompleteness, which we 
will consider in turn. 

First is the sampling completeness, meaning the fraction of all 
likely group members in the photometric sample for which we ob- 
tained a spectrum. This is primarily dependent on magnitude and 
position, with no significant colour dependence. For the group sam- 
ple, we consider all galaxies that have a photometric redshift con- 
sistent with the final group redshift, within twice the 68 per cent 
confidence limits of the photometric redshift. The targeting com- 
pleteness within .Rrms is computed on an individual group basis, in 
three r magnitude bins (22-23.25, 23.25-24, 24-24. 75). Note that 
our ap proach here is different from that described in Balo gh et alj 
d201 ll) . For the field sample, we simply measure the fraction of 
galaxies with photometric redshifts between 0.8 and 1.0 and found 
within our GMOS fields of view, with a spectrum. A linear function 
of this fraction and r magnitude is fitted. 

A second source of incompleteness is the failure to ob tain a 
reliable redshift from a spectrum. As in iBalogh et al.l d201 ll) . our 
success rate is very high, > 80 per cent for i < 24. There is no 
significant colour dependence of this fraction. Many of the failures 
are likely z > 1.4 galaxies (where [Oil] is redshifted out of the 
observed wavelength range) and thus outside the redshift range of 
interest here, so we apply no correction for this. 

Although the photometric redshift pre-selection greatly im- 
proves the efficiency of our observations, it introduces a potential 
bias. In particular, the process is sensitive to catastrophic failures, 



where objects in which the photometric redshift differs from the 
true redshift by many times the estimated uncertainty. The contri- 
bution of these outliers will not be appropriately recovered with 
our completeness weighting. To address this, in Figure Q] we show 
the difference between the spectroscopic and photometric redshift 
for all confirmed group members (see § 13. It with secure spectro- 
scopic redshifts and 0.7 < z p hot < 1.5, as a function of their 
(V ~ z) colour. Red points identify "priority 2" galaxies, which 
are those lower priority targets with z p hot more than 2<r away from 
the group redshift, but still within 0.7 < 2 p hot < 1.5. These seven 
galaxies contribute ~ 4.0 ± 1.5 per cent to the final group sample. 
From the present, small sample, these galaxies have colours repre- 
sentative of the rest of the group population. Moreover, only two 
of them (representing 1.2 per cent of the group sample) have stel- 
lar masses greater than the limits to which we consider our results 
complete (see § 13.4. 2t . Generally, however, we sample the priority 
= 2 population more sparsely than our primary sample. The sam- 
pling completeness of priority = 2 galaxies is ~ 0.15, or about 
0.2 times that typical of the top priority population, with a range 
of 10-30 per cent, depending on the group. Thus, we might expect 
that even after completeness corrections described above have been 
applied, we may be missing ~ 6 per cent of the group population 
above our completeness limit, due to this bias. 

Finally, when computing volume-normalized quantities like 
the stellar mass function for the field sample, we include a weight 
to account for the fact that the lowest luminosity objects are drawn 
from a smaller volume. 
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l°g,o(M lim /M ) 

Figure 3. Each curve shows the cumulative distribution of limiting stellar 
mass (i.e. at a given stellar mass, the fraction of the M/L r distribution that 
will be included in the sample) for a population of fixed (V — z) 0,9 colour 
and redshift. This is derived from the luminosity limit and the distribution 
of M/L r (see Fig[2} for that population. Weighting galaxies of that mass 
by the inverse of this number then accounts statistically for galaxies with 
high M/L r ratios that may be absent from the sample. The blue curves 
represent blue galaxies (V — z) 0,9 ~ 1.7, while the red lines correspond 
to (V — z) 0,9 ~ 2.9. Solid and dashed lines represent the lower and upper 
redshift range of the data, respectively (z = 0.8 and z = 0.95). 



3.4.2 Mass Completeness 

The r-selection of our survey introduces strong colour-dependence 
in the stellar mass completeness. In Figure [2] we show the stellar 
mass-to-light ratio as a function of (V — z) 0,9 colour, where the 
luminosity L r is computed from the selection magnitude (r) and 
the luminosity distance, without any k-correction. Therefore, the 
M I ' L r ratio directly relates the selection magnitude to the physical 
quantity of interest here (stellar mass), with distance as the only 
other relevant parameter. We divide the galaxies into those at the 
lower end of our redshift range (0.8 < z < 0.85, black points) and 
at the upper end (0.95 < z < 1.0, red points). 

The strong correlation with colour primarily reflects the dif- 
ferent star formation histories. At a fixed colour, there is a range 
of M/L r ratios that reflects a variety of effects, including physical 
effects like dust and metallicity variations. There are also varia- 
tions in the rest-wavelength coverage, which depends on redshift 
and SED shape. Finally, uncertainties in the SED-fitting procedure, 
resulting from photometric uncertainties, or the random sampling 
of bursts in the population modeling for example, will also con- 
tribute. The lines show the median M/L r at the lower (solid line) 
and upper (dotted line) ends of the redshift range. There is no sig- 
nificant difference between the two. The luminosity limit of our 
sample is well defined, by the r < 24.75 sample selection. This 
limit then corresponds to a range of stellar masses that depends 
strongly on colour. At fixed (V — jz) 0,9 colour, the maximum value 
of M / L r determines the 100 per cent mass completeness limit. For 
example, the reddest galaxies have M/L r as high as 30, and thus 
the sample is only fully complete for M sta r > 1O 1O,8 M0. 

However, this is unnecessarily conservative. With the assump- 



tion that the M/L r distribution, at fixed (V — z) ' 9 colour, is not a 
strong function of stellar mass, we can use the distribution shown in 
Figure|2]to statistically correct for the missing, high M/L r galax- 
ies. We construct the cumulative distribution of limiting masses, 
given by the r-band luminosity limit at a given redshift multiplied 
by the Mj 'L r , in bins of colour and redshift. An example is shown 
in Figure [3] Each curve shows, at a given mass, what fraction of 
the M/L r distribution (for galaxies of that type and redshift), will 
be included in the sample. Weighting galaxies of that mass by the 
inverse of this number then accounts statistically for galaxies with 
high Al/L r ratios that may be absent from the sample. If we are 
willing to accept statistical corrections as large as a factor ~ 5, we 
can extend the mass limit of blue galaxies by a factor ~ 4, and of 
red galaxies by a factor ~ 2. 

This correction is not strictly a statistical weight, but a sys- 
tematic correction. It means we are assuming that the distribution of 
dust, metallicity, and SED-fitting uncertainties is mass-independent 
for galaxies of a given colour over a relatively small mass range. 
This does mean that the galaxies we are missing are preferentially 
those with high M/L r ratios, e.g. those with more dust than aver- 
age. These may not be representative of the whole galaxy popula- 
tion at that mass and will remain a caveat to our methodology. 

3.5 Galaxy Sub-Populations 

It is now well known that most galaxies tend to populate either the 
"passive-sequence" of galaxies with negligible present star forma- 
tion, or the "blue cloud" of star-forming galaxies, with star forma- 
tion rates SFR oc M s " ar , where a is generally negative but n ot too 
far from unity (e.g. lNoeske et al.l2007tlwhitaker et al.l2012bl) . It is 
of interest to consider our sample in terms of these two populations 
and any galaxies that might be intermediate between the two. 

While this choice could be made based directly on estimates 
of the SFR, there are advantages to considering a classification 
based on more directly observ able quantities. It has been shown 
by us (iBalogh et alj|200^.l201 ll) and oth ers (e.g. lLabbe et alj2005l : 
IWolf et alJl2005l ; IWhitaker et alj|2012al) that the active and passive 
populations are well-separated in a colour-colour plane, as shown 
here in Figure [4] The x-axis colour, (V - z)°' 9 , brackets the 4000 
A break and is a good indicator of luminosity-weighted age, while 
the y-axis colour (J — [3.6] ) 0,9 , is sensitive to dust extinction. We 
plot all galaxies with good redshifts 0.8 < z < 1 within a broad 
mass range il/ star > 1O 9 ' 8 M0. The figure includes all group and 
field galaxies within the GMOS fields-of-view. We highlight with 
purple circles galaxies d etected in deep MIP S data, with 24/im flux 
greater than 0.071 mjy dSanders e t al. 2007). 

Figure[4]demonstrates the well-known fact that selecting pas- 
sive galaxies based only on an optical colour straddling the 4000 
A break (e.g. (V — z) ' 9 > 2.4) would include a significant pop- 
ulation of dust-reddened galaxies with significant star formation. 
Therefore, using optical data alone in the colour-colour plane, we 
seek to identify truly passive galaxies with no contamination from 
very dusty galaxies. First, a line is fit to the main distribution of 
galaxies, shown in Figure[4]as the blue dotted linfl The red dotted 
line in Figure|4]has an identical slope, but with an offset of —1.02 
mag to fit the second concentration of galaxies observed. Next, the 
distance of each galaxy from the two main distributions is mea- 
sured as the perpendicular distance from these lines. The distribu- 

3 The equation of this line is y = 0.89x — 0.28, where x = (V — z) 0,9 
and?/ = (J — [3.6]) - 9 . 
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Figure 4. Galaxies at 0.8 < z < 1 are shown in the colour-colour plane, 
aimed at separa ting dust-reddened spiral galaxies from the truly passive 
population (e.g. Balogh et a l. 2009, 2011). Confirmed group members are 
shown as filled symbols, while field members within our GMOS fields-of- 
view are shown with open symbols. The fits to the red and blue populations, 
are shown in the dotted lines. The thick green lines represent the boundaries 
defining the intermediate, green population. Galaxies with 24 micron MIPS 
detection within 3" are outlined in purple. 



tion of these distances is shown in Figure[5] with a clear separation 
between the two peaks evident. We also observe an intermediate 
population, located in between the mai n distributions. We w ill re- 
fer to these as the "green" galaxies, as in Balogh et al.l d201 ll) . After 
fitting the three populations with Gaussian functions, we define the 
boundaries of the green population so it has the minimal overlap 
with the passive and star-forming sequence. This is shown by the 
dotted lines in Figure|5] The resulting green box in Figure ^repre- 
sents our definition of the "transition" population. Note that there 
are additional vertical limits drawn using in the (V — z) ' 9 axis, in 
order to prevent very blue or very red objects from being identified 
as green galaxies. Si milar attempts to id entify green g alaxies have 
been m ade recently bv lGriitzbauch et al1 ( l201 lh and Whitake r et alj 
though the definitions differ from ours. 
In Figure [6] we show the (V — z) ' 9 colour dependence on 
stellar mass. Red, green and blue galaxies as defined above are rep- 
resented with like-coloured symbols. The straight lines represent 
the adopted stellar mass limit, defined as the mass corresponding 
to our r — 24.75 mag selection limit at the M/L r ratio for which 
50 per cent of our galaxies at that colour are expected to have lower 
M/L r . The solid and dotted lines represent the limits at z = 0.8 
and 2 = 1, respectively. We also show, as crosses, the seven group 
galaxies that originated from the priority = 2 sample. Most of these 
are low-mass galaxies, below the limit where we are complete ex- 
cept for the very bluest systems. 



3.6 Star Formation Rates 

While SFR can be estimated from our SED-fits, these can be dif- 
ficult to interpret, especially for low-SFR galaxies which are un- 
detected in the rest-UV. The SFR can instead be obtained directly 



c 




3 
Z 



0.5 1 

Colour offset from blue sequence 

Figure 5. The distribution of group galaxies positions (with masses 9.8 < 
log(Mstar/Mo) < 11.5) in the colour-colour plane, measured as the per- 
pendicular offset from the peak of the blue sequence. The histogram data 
includes the relevant statistical weights. Also plotted are Gaussian fits to 
the red, blue, and green population, and the summed function of all three 
fits. The vertical lines represent our identification of the intermediate, green 
population. 



from the [Oil] feature, as its luminosity is coupled to the H II re- 
gions, the sites of star-formation inside galaxies. However, this line 
is more strongly affected by dust extinction and metallicity than 
comparable indicators, like the Ha recombination feature. To com- 
pensate, a recent paper by Gil bank et al.l |3o 1 1 presents an em- 
pirical conversion between observed [On] luminosities and SFR at 
z — 0.1 in a statistical manner, calibrating the result to values de- 
rived using Ha and UV data: 



SFR 



J [o II] 



(4) 



M yr- 1 3.80 x 10 40 erg s" 1 a tanh[(a; -b)/c]+d' 

with a = -1.424, b = 9.827, c = 0.572, d = 1.700, x = 
log(M star /M ). 

We can measure the [Oil] luminosity by multiplying the 
equivalent width measurement with the luminosity of the contin- 
uum at its redshifted position. In the range 0.8 < z < 1.0, the 
r and i bands bracket the [Oil] line and the continuum luminos- 
ity is estimated by interpolating between the total magnitudes in 
those bands, measured from the broad-band images. As long as the 
locally-calibrated relation in Equation [4]holds at z = 1, these SFR 
are dust-corrected. However, this only includes the average dust 
correction for galaxies of that stellar mass and does not include any 
individual correction based on the observed properties. 

A comparison of the [On] - derived SFR using this method 
and the FUV + IR SFR for our sample show a systema tic normal- 
ization offset between th e two, as noted by others (e.g. IPatel et al] 
|20T it iKarim et al J 1 20 1 II) . It is not clear which method gives the 
right answer, but in principle the FUV + IR should be the most di- 
rect, with the IR bolometric correction being the largest systematic 
uncertainty. At this redshift range, 0.8 < z < 1.0, the IR lumi- 
nosity estimated fro m the 24 micron line should remain accurate 
dMagdisetal.ll201ll). 

The iGilbank et alj fcoich calibration would be expected to 



© 0000 RAS, MNRAS 000, 000-000 



Satellite quenching in z ~ 1 Galaxy Groups 9 



Log Stellar Mass 


Blue Group 


Blue Field 


Green Group 


Green Field 


Red Group 


Red Field 


9.5 - 9.8 


21 


26 










9.8- 10.1 


17 


50 










10.1 - 10.4 


20 


28 


2 


2 






10.4- 10.7 


12 


16 


7 


1 


14 


4 


10.7- 11.0 


13 


15 


4 


5 


11 


6 


11.0- 11.3 


1 


5 


2 


2 


11 


7 


11.3- 11.6 










1 


1 



Table 3. The number of galaxies in each mass bin, based on the colour definitions in Figure [4] The count is 
restricted to galaxies in the GMOS fields-of-view and with total weights (spectroscopic and mass completeness) 
of between and 5. Within the GMOS field of view, for 73 of the 304 galaxies in the sample, the zCOSMOS 
spectra were used. 
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Figure 6. The (V — z) 0,9 colour is shown as a function of stellar mass, for group galaxies (left panel) and for the field popluations (right panel). Galaxies 
are divided into three types, based on their location in colour-colour space (Figure |4}- Red, green and blue galaxies are colour-coded appropriately. The solid 
line and dotted line represent the stellar mass limit represented by the median M/L r ratio and our r = 24.75 magnitude selection, at z = 0.8 and 2 = 1, 
respectively. Galaxies with 24 micron MIPS detection within 3" are outlined in purple. The few cross symbols in the left panel correspond to group galaxies 
from our more sparsely sampled, priority 2 population, in which the z p hot is more than 2cr away from the group redshift. Most of these are low-mass (faint) 
galaxies, where even the primary sample is highly incomplete, except for the bluest galaxies. 

underestimate the true SFR if these relatively massive galaxies at 
z = 1 are dustier than galaxies of similar mass at z = 0, which 
may not be unreasonable. Recent studi es have found that th e aver- 
age dust attenuation may peak at z ~ l jCucciatietalj2012h . Thus, 
we elect to multiply the [Oil] - derived SFR by a factor of 3. 1 in all 
subsequent analysis, as explained in Appendix [A] The difference 
between these two methods will be explored in a subsequent pa- 
per, Mok et al. (in prep), which will also compare the results from 
other star formation indicators on our sample, including the use of 
SED-fitting. 

4 RESULTS 

4.1 The Correlation between sSFR and Stellar Mass 

In Figure [7] we show the specific star formation rate as a function 
of stellar mass, considering galaxies within the GMOS fields of 
view and with secure redshifts 0.8 < z < 1. Galaxies are colour- 



coded according to their classification into red, green, and blue, as 
described in § 13.51 The solid and dotted lines represent the com- 
pleteness limits at z — 0.8 and z = 1, respectively. 

We observe that the blue galaxies form a well -defined se- 
quenc e, in good agreement wit h previous results from lElbaz et ail 
d2007h at the high mass end and Whita keret alj J2012ah above their 
completeness limit of Af s tar ~ 1O 1O M0. Our data extend this trend 
to lower masses, M ata r ~ 10 9 ' 5 Mq. In addition, we observe that 
the green galaxies are not well separated from the red population, 
as the [Oil] feature is rarely detected in either population. A small 
number of blue galaxies also do not have a detected [Oil] emission. 
At the high mass end, this may be caused by the [Oil] luminos- 
ity being dominated by the continuum component and difficulties 
in the detection of relatively weaker [Oil] feature from the spec- 
tra. Variations in dust extinction from the average, and the large 
uncertainties especially for fainter galaxies, will also contribute to 
variations from a simple correlation between colour-class and line 
strength. 
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Figure 7. The specific star formation rate measured from the [On] feature is 
shown as a function of stellar mass for galaxies at 0.8 < z < 1 with secure 
redshifts. Group members are shown as filled symbols, and field members 
within our GMOS fields-of-view are shown as open symbols. Only galax- 
ies with total weight (spectroscopic and mass-completeness) of between 1 
and 5 are included. The solid line and dotted lines represent our 50 per cent 
completeness limit due to the M/L r distribution at fixed sSFR, at z = 0.8 
and 2 = 1, respectively. Points are colour-coded according to their clas- 
sification as red, blue, and green. Galaxies with significant X-ray sources 
within 5" are outlined in black, while galaxies with 24 micron MIPS de- 
tection within 3" are outlined in purple. The star-forming sequence from 
lElbaz et al 1 l2007n ( thick pink dashed line ) and Whitaker et al. ( 20 1 2j|) ( thin 
pink dashed line) at z = 0.9 are included for comparison, up to their stated 
mass completeness limits. Note that those galaxies with a negative [On] 
equivalent width or with [Oil] detections less than the measurement error 
are converted into upper limits. 



The slope of the sSFR-stellar mass relationship potentially 
holds important clues about the processes regulating galaxy for- 
mat ion, and can put v ery strong constraints on feedback models 
(e.g lBower et al ] |2012l) . However, the measured slope will depend 
sensitively on the completeness limits and the definition of star- 
forming galaxies. In particular we are certainly missing galaxies 
with substantial SFR at M s tar < 1O 1O M0, which would mean 
the overall slope of the star-forming sequence is likely to be flat- 
ter than observed here, despite the inclusion of statistical weights 
to partially offset this. 

In Figure [8] we show the weighted median sSFR for blue and 
green galaxies, stacked in equal mass bins, with errors calculated 
using the bootstrap method. For the most massive blue group bin, 
which contains only one galaxy, the error of the measurement is 
shown. For the green galaxies, cases where the weighted median 
would indicate an undetected [On] emission are shown as upper 
limits. We observed that for M star > lO ia3 M , the blue and 
green field populations show slightly higher sSFR than their group 
counterparts. Below M s tar ~ 1O 1O ' 3 M0, the blue group and field 
population do not show a statistically significant difference, which 
may be caused by incompleteness in our sample, since we could 
only detect galaxies with high sSFRs in this mass range. 

To determine the cause for this possible diver gence in the 
sSFR of group and field population, we used the Iwhitaker et alj 
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Figure 8. The binned median sSFR is shown as a function of stellar mass for 
galaxies at 0.8 < z < 1, within the GMOS fields-of-view and total weights 
(spectroscopic and mass-completeness) between 1 and 5. Group galaxy bins 
are shown with filled symbols and field galaxies as open symbols. The solid 
line and dotted lines represent our 50 per cent completeness limit due to 
the M/L r distribution at fixed sSFR, at z = 0.8 a nd 2 = 1, respectively. 
The star-forming seque nce from lElbaz etafl (2007) (thick dashed line) and 
Whitaker et al. (2012a) (thin dashed line) at 2 = 0.9 are included for com- 
parison, up to their stated mass completeness limits. The blue points are 
the binned values for the blue galaxies and the green points indicate green 
galaxy bins. For the green galaxies, cases where the weighted median would 
indicate an undetected [Oil] emission are shown as upper limits. Note that 
for the most massive blue group bin, which contains only one galaxy, the 
error of the measurement is shown. 



d2012al) relationship between stellar mass and sSFR at z = 0.9 
and calculate the difference between their sequence and the indi- 
vidual galaxies' sSFR. The distribution of the residuals for the high 
mass blue galaxy population (A/ S tar > 1O 1O ' 3 M0) is shown in Fig- 
ure [9] including the relevant spectroscopic and mass-completeness 
weights. Qualitatively, the main difference between the two shapes 
is the lack of high sSFR galaxies in the group sample. Quantita- 
tively, the two distributions can be compared by calculating their 
weighted means, with the associated error provided using the boot- 
strap method. The field population has a mean log residual of 
0.27 ± 0.1, while the mean residual for the group population is 
— 1.48 ± 1.6. Performing the standard Kolmogorov-Smirnov (KS) 
test on the unweighted sample shows that the group and the field 
have a 17% chance of originating from the same distribution. We 
conclude, therefore, that the small offset we observe between the 
star-forming main sequence of z ~ 1 group and field galaxies is 
not statistically significant. 

4.1.1 Stellar Mass Functions 

In this section we consider the stellar mass function of all galax- 
ies in our spectroscopic group an d field samples. We will com- 
pare extensively with the results of lGiodinTe t al. (2012), who pub- 
lished mass functions for a similar (and partially overlapping) 
sample of groups, based on photometric redshifts only. Benefit- 
ing from their larger sample, they were able to further divide the 



© 0000 RAS, MNRAS 000, 000-000 



Satellite quenching in z ~ 1 Galaxy Groups 1 1 



i 1 1 1 1 i 

Shaded = Group 
Unshaded = Field 




•1 1 

log sSFRj „] Residuals 



Figure 9. A histogram showing the residuals to the Whitaker et al.l d201 2al) 
line at z = 0.9, shown in Figure [8] for the high mass blue galaxy pop- 
ulation (M s tar > 10 10,3 Mq). The data includes the spectroscopic and 
mass-completeness weights, drawn from the same population as Figure [7] 
The histogram for the group galaxies is shaded, while the field galaxies are 
unshaded. 
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Figure 10. The stellar mass function for all group galaxies in our sample is 
shown as the solid black circles, and with lower limits indicating where our 
sample is incomplete for red galaxies. Open symbols represent the field sam- 
ple within the GMOS fields of view, slightly offset horizontally for clarity. 
The filled open purple triangles represent massive (M s tar > exlO^M^ ) 
and low-mass groups from the photometric analysis of Giodi nTet all i2012l) . 
while the dashed line represents their field sample. The Giodini data are all 



arbitrarily renormalized to match our data at M s t a 



10 1 
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group sample into low mass and high mass subsamples (divided at 
Mstar = 6 x 1O 13 M0). As groups in our sample come from both 
sides of this division, we compare our data to both results. In addi- 
tion to the deep spectroscopy, the novelty of our present work is the 
colour division of the galaxy population into red, blue, and green. 
By contrast, Giodini et al. had only divided their sample into two 
populations based on SED-type: passive and star forming. 

We begin by considering the stellar mass function of all group 
galaxies, computed by counting the number of galaxies in a given 
stellar mass bin, and applying the weights discussed in § l3.4l to cor- 
rect for sampling and M/L r incompleteness. We normalize each 
bin by the number of groups in our sample that potentially con- 
tribute (depending on their redshift), so the mass function shows the 
average number of galaxies per group. Then, we restrict the mass 
range so that all bins have contribution from at least four groups. 

In general, only the lowest-mass bin includes contributions 
from fewer than half our groups. The uncertainties in this and all 
subsequent mass functions are dominated by the uncertainties in 
the stellar masses themselves, which are systematic in nature. We 
account for these by simply comparing the mass function computed 
with all stellar masses at their la upper and lower values. 

The resulting mass function is shown in Figure[l0] as the solid 
black circles with error bars. The two lower limits correspond to 
measurements in mass bins for which our sample is incomplete 
for red galaxies, even after weighting. The open circles show the 
corresponding mass function for the field population, considering 
only the sample within the GMOS fields of view. These are ar- 
bitrarily renormalized to match our "per-group" normalizations at 
Mgtar ~ 10 10 ' 8 Mq. The filled and open purple tri angles represent 
massive and low-mass groups from lGiodini et al.l ( 120120 . while the 
dashed line represents their field sample. Again, all are renormal- 
ized to match at M ~ 10 10 ' 8 Mq . We will retain this normalization 
in subsequent plots, where the galaxy population is divided by type. 
The mass function shapes are generally in excellent agreement over 
the full mass range, though the mass funct ion of our fie ld sample 
appears to be somewhat flatter than that of Giodini et^al] d2012l) . 



4.1.2 Type-Dependent Mass Functions 

We now divide the stellar mass functions by galaxy type as defined 
in § 13.51 and start with the red population in Figure QT| We only 
show data to the limit where our completeness drops below 20 per 
cent. Above that, the data is weighted to correct for incompleteness. 

The passive galaxies in our groups exhibit a steep mass func- 
tion, with no sign of a turnover above our limiting stellar mass of 
-Mstar = 10 10 3 Mp). Apart from the lowest-mass bin, our data is in 
excell ent agreement with the photometric analysis of Giodi ni et alj 
J2012I) , with the relative normalization based on the full sample 
shown in Figure [10] This difference in the l owest-mass bin is also 
observed between the mass function in the iKnobel et alj d2012bh 
sample, which is based on an average of a large number of sparsely- 
sampled groups with a wide range of masses. Our work is comple- 
mentary because we have a robust sample of groups with reason- 
ably well-determined masses, but our statistics are not yet good 
enough to make a strong claim about this lowest mass point. 

In contrast with the groups, our field sample shows a lack of 
passive galaxies below Af s tar = 10 10 ' 7 Mq. This is reflecting the 
difference apparent in Figure [6] where the low-mass end of the red 
sequence of the field sample is apparently less well-populated than 
the group sample. This turnover is more prominent than in Giodini 
et al.'s field sample. The difference may be a cosmic variance ef- 
fect, as the area of our survey is smaller, or it could reflect the fact 
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Figure 11. The stellar mass function for passive group galaxies in our sam- 
ple is shown as the filled red points; the field sample is shown with open 
circles. This is again compared with the high- and low-mass photometric 
group samples lGiodini et al.M2012l) . shown as the filled and open triangles, 
respectively, as well as their field (dashed line). The normalization of the 
Giodini et al. data are all fixed according to the total mass functions. We 
find a steep stellar mass function for passive galaxies, with no indication 
for a turnover in the group population. 
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Figure 13. The fraction of passive galaxies in both group and field samples 
is shown as a funct ion of stellar mass. W e compare these results to the pho- 
tometric sample of Gio dini et al 1l l2012l) . shown as the purple triangles with 
the same meaning as previous figures. The blue lines indicate the observed 
relations in the loc al universe as measu red from the SDSS using a very sim- 
ilar method, from lMcGee et al] l l201 lh . The solid line represents the local 
group population, while the dashed line is the field. 
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Figure 12. As Figure[TT] but for the star-forming galax y population, s hown 
in blue. Purple triangles represent the groups from lGiodini et alj J2012h as in 
previous figures; note that their population of star-forming galaxies includes 
our green galaxies, which may explain the discrepancy at higher masses. 
The Giodini et al. field mass function for blue galaxies is shown as the 
dashed line. 



that our field definition is systematically underdense because the 
dominant group in each pointing has been removed. 

Figure [T2] shows the stellar mass function for the blue (star- 
forming gal axies) in our group sample, again compared with the 
field and the lGiodini etailhoiA data. The shape of the group mass 
function is very similar to that of the field, but with a normalization 
that is lower by a factor ~ 2. As these mass functions are normal- 
ized so the total mass functions match, this reflects a lower fraction 
of blue galaxies in our groups, at all stellar masses. 



4.1.3 The Fraction of Passive Galaxies 

The steep mass function for passive galaxies leads to a large frac- 
tion of such galaxies at low masses. We show this in Figure [13] 
by dividing the passive-galaxy mass function by that for the full 
sample. We compare this with the fraction in our spectroscopic 
field sample ( open symbols), and the photometric field sample of 
lGiodinietai] ( l2012h . At high masses, the two are indistinguishable, 
though the uncertainties are large given our small volume. But at 
A/star < 10 11 Mq the groups remain dominated by purely passive 
galaxies, at ~ 60 per cent of the population, while the field is dom- 
inated by star-forming galaxies. Again, this reflects the absence of 
low-mass passive galaxies apparent in Figure [6] 

Comparing ne xt with the photometric group sample of 
iGiodini et al] d2012l) . we see good agreement over most of the mass 
range. However, for AW < 1O 1O ' 6 M we find the group pop- 
ulation remains dominated by passive galaxies, where Giodini et 
al. note a continuing decline. The origin of this discrepancy is not 
clear. It possibly indicates a bias with the photometric method of 
Giodini et al., but as the most significant discrepancy is in a single 
bin, for which passive galaxies are only detectable in our lowest- 
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redshift groups, this result should be tested with deeper spectro- 
scopic data and a larger sample. 

As it stands, our data suggest that groups at z ~ 1 were al- 
ready dominated by passive g alaxie s at all stellar masses M star > 
10 10,3 M Q . In iMcGee et alj feOllh we presented an analysis of 
the SDSS and our GEEC groups at z ~ 0.5, and identified pas- 
sive galaxies based on their SED-derived SFRs. The results from 
that SDSS analysis are shown as the solid and dashed line for 
the group and field samples, respectively. Our field sample clearly 
differs from that at z ~ 0.1, with a steeper dependence of this 
fraction on stellar mass and an almost negligible passive popu- 
lation at the lowest masses. On the other hand, the passive frac- 
tion in the 2 = 1 groups appears to be highest for the lowest 
masses, where the observed fractio n is similar to that observed lo- 
cally (see also iGeorge et al.ll201 lb . Again, much depends on the 
bin with the lowest stellar mass galaxies, and a deeper, larger spec- 
troscopic sample is required to confirm this surprising result. Cer- 
tainly, we do not find any evidence that the fraction of passive group 
galaxie s at a given s t ellar m ass is ever lower than in the field. Us- 
ing the Knobel et al. ( 2009) catalogue of the spectroscopic groups, 
llovino et alj ( 120 100 have identified a mass-dependent environmen- 
tal effect on the fraction of blue galaxies up to z ~ 0.8, measured 
using their (U — B) rest frame colour. This effect is strongest in 
the stellar mass range M s t ar < 1O 1O ' 6 M0, which is also the range 
where we observe the largest difference in the fraction of passive 
galaxies between the group and fi eld. Using the local d ensity esti- 
mator for the zCOSMOS sample, ICucciati et alj d2010h have also 
found that the local density modulates galaxy colour for this in- 
termediate mass range (10 la2 M Q < M atar < 10 10 7 M Q ). Our 
spectroscopically complete sample demonstrates this density de- 
pendence was in place at z ~ 1, corresponding to a look back time 
up to 1 Gyr earlier than the maximum probed in these studies. 



4.2 Spectroscopic Properties of the Red, Blue, and Green 
Populations 

Next, it is instructive to consider the average properties of each 
population. Figure [14] shows an example of the stacked spectra for 
the six sub-populations, with relatively high signal-to-noise ratios 
of between 7.0 per A (in the range of 4050 A to 4100 A) for the 
green field stacked spectra to 12.7 per A for the more numerous 
blue group stacked spectra. 

The average spectra confirm that green galaxies have features 
that are intermediate between blue and red galaxies. They have 
less prominent absoiption lines, and do not resemble typical post- 
starburst type galaxies. Second, [Oil] emission is present at a low 
level, much more so in the field population. Thus, some star forma- 
tion is still ongoing in these galaxies. 

Using the stacked spectra, the "average" spectral properties for 
the red, blue, and green populations can be determined by measur- 
ing the [Oil] and US strengths. Note that our stacked spectra con- 
tain enough signal to distinguish the presence of weak emission- 
filling in the US absorption line of blue galaxies. We correct for 
this in the calculation of the average US equivalent width, by fitting 
multiple Gaussian functions. The equivalent width of the emission 
line is then added to the total equivalent width obtained using the 
automated bandpass method, as described in §[330 For the blue 
population, we separated the populations into a low mass and high 
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Figure 14. Rest-wavelength stacked spectra for galaxies at 0.8 < z < 1 
and 9.5 < log(M s tar/M©) < 11.6, including the spectroscopic and 
mass-completeness weights, for the six main classification used in this pa- 
per. The spectra are arbitrarily normalized within the wavelength range of 
4050 — 4100 A. The location of the important lines (in order: [Oil] , Cal- 
cium H&K, and H<5) are also indicated on this plot. The figure on the right 
show an expanded section of the spectra around the H5 line. 



mass sample at log(Af sta r/Mo) = 10.3. This decision was mo- 
tivated by the differences observed in the distribution of sSFR in 
§ 14. ll between the high mass and low mass blue galaxies, caused in 
part by the onset of incompleteness. 

First, the sSFR for these four bins are shown in Figure [15] 
with filled and open circles for the group and field populations, re- 
spectively. For the red and blue populations, this average shows at 
most a weak dependence on environment. For the blue galaxies in 
particular, where the statistical uncertainty on the stacked spectra is 
negligible, we observe a significant difference between the average 
sSFR of the low and high mass blue sample. While the group galax- 
ies have only ~ 3 per cent lower sSFR than the field sample for the 
low mass blue population, the difference is ~ 55 per cent for the 
high mass sample. The apparently high significance of this differ- 
ence is misleading, however, as the error bar only reflects the fact 
that the [Oil] feature is measured with high signal-to-noise in the 
stacked spectrum. Intrinsic variations in the line strength amongst 
galaxies means the average SFR is not really so well determined, 
as already seen in Figures [8]and [9] 

Interestingly, the intermediate, green galaxies show an even 
larger difference, as green group galaxies have an average sSFR a 
factor ~ 3 lower than the field, when measured using the stacked 
spectra. This result does need to be treated carefully, as this is a 
relatively small population defined by sharp boundaries amongst a 
continuous distribution (Figure [5}- Both statistical and systematic 
(i.e. the connection between colours and underlying star formation 
history) uncertainties mean that the green population likely con- 
tains "contamination" from the adjacent blue and red populations. 
Since the red fraction is a strong function of environment, this can 



4 Note that for the lower mass blue galaxies bin, the emission component 
was measured to be 1.5 A for the group galaxies and 1.3 A for the field 



galaxies. For the higher mass bin, the emission component was measured 
to be ~ 1.8 A for the group and field population. 
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Figure 15. The sSFR for the red, green, and blue populations. The galaxies 
are restricted to 0.8 < z < 1 and include galaxies within the GMOS 
fields-of-view and total weights (spectroscopic and mass-completeness) of 
between and 5. Bins of group galaxies are shown with filled symbols and 
field galaxies as open symbols. The red and green bins contain galaxies with 
masses 9.5 < log(Af s tar /Mq) < 11.6. For the blue population, another 
cut was made, separating the two populations at log(Af a tar / 'Af©) = 10.3. 
The circle symbols indicate the results from stacked-spectra measurement 
of [Oil] sSFR. The squares symbols are data from the same bin, except 
combined using the weighted median method. 

lead to an apparent environmental dependence amongst these green 
galaxies. 

When we compare their sSFRs from the weighted-median 
of the individual galaxy measurements (see § 13. 3t , shown as the 
squares in Figure [15] the difference between group and field still 
exists. However, the error bars now include the bootstrap contribu- 
tion to the error due to the scatter of the parent distribution. The 
larger error bars indicate that these differences in the sSFR are not 
statistically significant, given the small sample sizes and the wide 
spread in the properties of the binned galaxies. Therefore, we can- 
not confidently claim to detect any significant suppression of sSFR 
amongst group galaxies of a given spectral type, relative to the field. 

Next, we consider the H<5 absorption line. This spectral fea- 
ture is most prominent in A-type stars and is a signature of re- 
cent star formation, especially for galaxies with low instanta- 
neous sSFR. An especially strong H<5 signature, with an equivalent 
width greater than 5 A, may indicate a truncation in star forma- 
tion, due to the domination of intermediate age A-type stars (e .g. 
ICouch&Sharpleslll987l : iBalogh et alJI 19991: ILemaux et alj2012h . 

In Figure [16] we show the US equivalent width measured us- 
ing the stacked spectra as the circular points. One potential source 
of concern is that the corrections to the telluric absorption at the 
A-band may not be perfect. They often leave significant residu- 
als which unfortunately overlap with US over the redshift range, 
0.818 < z < 0.895. To ensure that this contamination is not 
affecting our results, we remove galaxies in the redshift range 
0.818 < z < 0.895 and show our updated results for this restricted 
sample with triangular points. 

We find differing environmental signatures in US between the 



Figure 16. The H<5 equivalent widths for the red, green, and blue popula- 
tions. The galaxies are restricted to 0.8 < z < 1 and include galaxies 
within the GMOS fields-of-view and total weights (spectroscopic and mass- 
completeness) of between and 5. Bins of group galaxies are shown with 
filled symbols and field galaxies as open symbols. The red and green bins 
contain galaxies with masses 9.5 < log(M s tar/M0) < 11.6. For the 
blue population, another cut was made, separating the two populations at 
log(M s tar/A40) = 10.3. In addition, the measurement of the blue galax- 
ies includes a correction for H<5 emission, computed using a Gaussian fit to 
the stacked spectra. The triangles represents the removal of galaxies in the 
redshift range 0.818 < z < 0.895, in order to avoid the telluric-dominated 
spectra. 

red, green, and blue populations, though they are not highly sig- 
nificant given our uncertainties. The blue galaxies (M star / Mq > 
10 10 ' 3 ) show a weaker absoiption for the group sample, as com- 
pared to the field, which may be related to the observed differences 
in their sSFR from Figure [15] However, this significant difference 
disappears when we consider the sample unaffected by telluric ab- 
sorption, so we again treat this result with caution. Another result to 
note from the US measurements is that the green population is not 
dominated by post-starburst galaxies (VKc^HcJ )> 5 A), as observed 
from their spectra. Surprisingly, the H5 strength of the green popu- 
lation is much less than that of the blue population, which suggests 
that the last period of significant star formation occurred ^1 — 2 
Gyr ago. In addition, the green field galaxies have a weaker H<5 line 
than even the red galaxies, though this result should be treated with 
caution. The systematic errors are not included in the error bars, as 
demonstrated by the substantial shift after the removal of the tel- 
luric dominated galaxies. Furthermore, these green galaxies may 
have an undetected US emission component in their low signal-to- 
noise spectrum, especially for the field population in which we see 
significant [Oil] emission. 



5 DISCUSSION 

We proposed in Balo gh et alj d201 lh that the intermediate green 
population may re present a transition from the blue cloud to re d 
sequence (see also iGeorgakakis et al]|2008l : ILemaux et alj|2012l) . 
In Figure [17] we show the ratio of the "green" and "blue" mass 
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Figure 17. The ratio of the mass functions of the transition population and 
the star-forming population are shown, for our group (filled circles) and field 
samples (open circles). The ratio increases with stellar mass, and shows 
little environmental dependence. 



functions. If the former population is indeed a transition from the 
latter, this ratio should give some indication of the rate at which 
this transformation is occurring. What we see is that this fraction 
is largest for the highest mass galaxies, and is ~ 20 per cent at 
Mstar < 1O 1O ' 7 M0. There is no strong evidence that the frac- 
tion is any higher within groups, though the uncertainties are large 
enough to be consistent with an interestingly-large difference. If the 
green galaxies do represent a transition population, many of them 
may be more closely a ssociated with the mass-quenching popula- 
tion jPeng etal.|[2O10l) . Another possibility is that some of them 
represent predominantly passive galaxies that have been partially 
rejuvenated, for example following a minor merger event which 
provides fuel for additional star formation. On the other hand, many 
of the field, green galaxies may be satellites of smaller groups or 
large galaxies. Thus the possibility remains that this population is 
largely environment-driven, a result that can be tested with highly 
complete spectroscopic surveys of the field at this redshift. 

Despite the fact that most of the transition population do not 
appear to be related to the group environment, the fact that we ob- 
serve significantly higher fractions of passive galaxies in the groups 
than the field, for A/ S tar < 10 10,7 Mq, implies that environment 
does play a role. Assuming that the properties of the observed field 
galaxies are similar to the pre-accreted population, we can com- 
pute an efficiency factor that describes the fraction of field galaxies 
at a given mass that would have to be quenched to ma tch the red 
fractio n at a given over-density. Following the spirit of IPeng et alj 
we define this environment quenching efficiency at stellar 
mass M s tar, as: 

_ /rod (group, M s tar) ~ /rcd(field, M sta r) ,„ 
/blue (field, Mstar) ' 

Sim ilar definition s are a dvocated by Ivan den Bosch et all 1 20081) 
and Wetzel et al.l J2012h . Using the results from McGee et all 
J20T l]) we compute e p at z = and z = 0.5, and compare this with 
our new data, in Figure Q~8] While the uncertainties are still fairly 
large, a value of e p ~ 0.4 is consistent with our data and is gener- 
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Figure 18 . The environment quenching efficiency, e p , as defined by 
IPeng et al] J2010h . is shown for four different stellar mass bins and three 
redshifts. The z = data are based on the SDSS; these and the z = 0.5 
data are taken from the measurements in McGee etalJlEoi 1 ), interpolating 
to the mass bins used here. 
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Figure 19. From the galaxy formation model of Bower et al. (20061). we 
show the cumulative fraction of satellites in groups of the indicated masses 
and at the indicated redshifts (2 = 0.0, z = 0.83, z = 0.98), as a function 
of the time since they were last identified as a "central" galaxy. The solid 
line indicates accretions into groups with log(M g r P /MQ) > f3.0 and 
dashed line indicates groups with log(Af g r P /AfQ) > f3.5. The model 
gives an indication of how long satellite galaxies may have been sensitive 
to environment-quenching. The stellar mass cut-off for the satellites is set 
at> 10 10 M Q . 



ally independent of both stellar mass and redshift. The value for the 
environmental quenching efficiency is also consistent with similar 
analysis done up to z ~ 0.8 from the larger zCOSMOS spectro- 
scopic group catalogue dKnobel et alj|2012bh . There are hints of a 
smaller e p for the higher masses, though the uncertainties are large. 



IWetzel etal.N2012l) recently argued, based on the shape of the 
SFR distribution amongst satellite galaxies in the SDSS, that the 
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shutdown of star formation in satellites is best characterized by 
a long (~ 3 Gyr), mass-dependent delay where evolution contin- 
ues as normal, followed by rapid (r q < 1 Gyr) quenching. If this 
timescale does not evolve with redshift, we might therefore expect 
that star-forming galaxies accreted between 3 Gyr and 3 + r q Gyr 
ago are in the process of quenching at the epoch of observation. 

However, the long delay, also advocated by |Pe Lucia et alj 

may be difficult to reconcile with the large environmen- 
tal effe ct we observe. The galaxy formation models of Bow er et"al] 
d2006h predict the fraction of satellites accreted for groups of our 
mass and redshift as a function of time, as shown in Figure [19] The 
results show that the fraction accreted in time At is ~ 20 At/ Gyr 
per cent of the total observed population. Less than ~ 30 per cent 
of satellites are predicted to have been a satellite for as long as 3.5 
Gyr before the epoch of observation. Thus it is difficult to see how 
an environmental signature as strong as that seen in Figure [T3lcould 
arise. T he longer timescale of ~ 7 Gyr advocated bv lDe Lucia et alj 
d2012l) is even more problematic, as this is longer than the age 
of the Universe at z = 1. If the "delayed-then-rapid" scenario of 
IWetzel et alj l2012h is correct, then it seems the delay must be con- 
siderably shorter at z = 1. For e xample, if it scales with the dynam- 
ical time as r oc (1 + z) _3/2 ^Tinker & WetzefeoiOl) . we might 
expect the relevant timescale at z = 1 to be 1-2 Gyr. This would 
likely be consistent with our observations. 

In principle we can use the observed fraction of transition 
galaxies to try and constrain the quenching time, r q . We find that 
green galaxies make up 10 per cent of the total galaxy population, 
at M s tar < 1O 1O ' 7 M0. This must be an upper limit to the fraction 
undergoing transition at the observed epoch, as at least some are 
likely to result from mass-quenching, or rejuvenated star formation 
in old galaxies. Given the high accretion rates shown in Figure [T9l 
the quenching time must be quite short, r q ~ 0.5 Gyr or less, in 
order to avoid overproducing galaxies in this phase of evolution. 
In a future paper, we will use a simple model, using Figure [19] to 
predict the red fractions and the expected H<5 strengths of the blue, 
green, and red populations using population synthesis models. This 
will allow us to constrain r q and the delay timescale, using our ob- 
servations. 



6 CONCLUSIONS 

We have presented highly complete GMOS-S spectroscopy for 
galaxies with r < 24.75 in the fields of 11 galaxy groups at 
0.8 < z < 1. We use an optical-NIR colour-colour diagnostic dia- 
gram to divide the galaxies into three subpopulations: star-forming 
(blue), passive (red), and intermediate (green). Using SED-fitting 
techniques and spectral analysis, we present a comparison of group 
populations with the field, down to a stellar mass of M sta r > 
10 9,5 A'/ Q . We conclude the following: 

• The strongest environmental dependence is on the fraction of 
passive galaxies. For galaxies with 10 10 ' 3 < Mstar/M© < 10 11 , 
groups are dominated by passive galaxies, which only constitute 
~ 20 per cent of the field at these stellar masses. 

• For the normal star-forming galaxies, the average sSFR for the 
high mass (Mstar / 'Mq > 10 10 ' 3 ) group sample is ~ 55 per cent 
lower than the field population, but the difference is not statistically 
significant. For the low mass blue galaxies, where our sample is 
incomplete, we observe no difference in sSFR between the group 
and the field. 

• The green galaxy population represents approximately the same 
fraction of the star-forming population at all stellar masses, in both 



the group and the field. Neither the group or field spectra show very 
strong US absoiption, which indicates that their populations are not 
dominated by post-starburst galaxies. Green galaxies in the group 
exhibit a factor of ~ 3 lower sSFR than the field, calculated using 
their stacked spectra, suggesting that the origins of these transi- 
tional systems may be environment dependent. 

In the context of the simple model proposed by Pen g et ail 

our results are consistent with an environment-quenching 
efficiency of ~ 40 per cent (for z < 1). There are hints that this 
value may be smaller for galaxies with higher stellar masses. Fur- 
thermore, the abundance of possible "transition" galaxies requires 
a relatively short quenching timescale of r < 0.5 Gyr. Various 
authors have argued there must be a long timescale, on the order 
of ~ 3 — 7 Gyr, associated with satellite quen ching dMcGee et alj 
l2009MDe Lucia et alj201 AlWetzel et al.l2012l) . This could be asso- 
ciated with a delay, before star formation is truncated on a shorter 
timescale. However, such a long delay time seems incompatible 
with the large environmental signature we see at z = 1, where most 
galaxies have not spent that long as satellites. Eit her this timescale 
must evolve so that it is shorter at higher redshift (Tinker & Wetzel 
I2010D . or perhaps only a fraction of galaxies ever respond to their 
environment. Better constraints could be obtained by extending this 
analysis to lower-mass galaxies with low or negligible star forma- 
tion, which will require deep integrations of NIR- selected samples. 
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APPENDIX A: NORMALIZATION DIFFERENCE 
BETWEEN DUST-CORRECTED [On] SFR AND FUV+IR 
SFR 

To calculate the FUV+IR SFR, we first k-correct the UV and optical 
photometric data from the zCOSMOS catal ogue using KCORRECT 
IDL software of Blanton & Roweis (2007) to their rest frame. To 
do this, we apply corrections to the filter that most closely matches 
the rest wavelength of the desired filter, rather than taking the 
synthetic magnitude from the template fit. This allows us obtain 
reasonable upper limits for galaxies that are undetected in rest- 
frame FUV flux. Next, MIPS 24 micron da ta is converted into to - 
tal IR luminosity with the templates from Chary & E lbazl d200lh . 
We then combine FUV and IR luminos ities into SFR estimates 
with the prescription of lHao et alj (201 1). Finally, galaxies with- 
out MIPS 24 micron data had dust attenuation estimat ed using the 
(FUV - NUV)° relationship from lSalim etail d2009l) . 

In Figure IA1I the log difference in the SFR between the two 
methods, logio (SFR[ FUV +ir]) - logio (SFR[ nj), is shown as a 
function of stellar mass. A comparison of the [Oil] dust corrected 
SFR and the FUV + IR S FR for our sample show a normal iza- 
tion offset, as noted in (e.g. lPatel et alj201 lHKarim et al]|201 lh . A 
mass-independent difference is observed, and so the [Oil] - derived 
SFR are multiplied by a factor of 3. 1 for the purpose of this paper. 
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Figure Al. Logio (SFR[ FUV+IH] ) - logio (SFR[ OII ] ) is plotted as a func- 
tion of stellar mass for galaxies at 0.8 < z < 1 in the GMOS fields- 
of-view. Group members are shown as filled symbols, and field members 
within our GMOS fields-of-view are shown as open symbols. Only galaxies 
with total weight (spectroscopic and mass-completeness) of between 1 and 
5 are included. Points are colour-coded according to their classification as 
red, blue, and green. Galaxies with 24 micron MIPS detection within 3" 
are outlined in purple. The black line indicates the mass-independent offset 
between the FUV + IR and Oil star formation rates. Note that those galaxies 
with a negative [Oil] equivalent width or with [Oil] detections less than the 
measurement error are converted into lower limits, while undetected FUV 
galaxies are shown as upper limits. 

This will be explored in more detail in a future paper (Mok et al., 
in prep.). 
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